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Outline 


Utilisation of micro-architectured structure as a heat sink medium. 

1. Analytical approach using two-equation model. 

2. Experimental results of pressure drop and heat transfer measurements 

Novel lightweight metal foams 

1. Examples employing metal foams as a heat exchanger 

2. Analytical approach using two-equation model. 

3. Experimental results of pressure drop and heat transfer measurements 
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Heat sink model using the LFMs 

Lattice-frame material (LFM)s are triangulated 3-D structures 
(Tetrahedral cell based) 

Greater lateral stiffness 
and lower cost per unit weight over 
honeycombs and metal foams 

- Surface area density : 123.68 [1/m] 

- Relative density ( or porosity) : 0.062 (0.938) Unit ceil of the lfms 

Forced air convection heat transfer measurements were conducted 
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Configuration of Lattice-Frame Materials 


Two orientations denoted as below were tested for heat transfer 
and pressure drop measurements 

O-B 



(b) Orientation-B (frontal view) (c) Crossflow arrangement of LFMs 
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Honeycomb 


0.127m 
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Experimental conditions & set-up 


Specification of LFM 

Typical operating parameter 

LFM cell bar diameter (d) 

0.002 m 

Inlet coolant mean velocity 

1.0-26.0 m/s 

Longitudinal cell pitch ( S x0 ) 

0.0127 m 

Bar Reynolds number 

120<Re, <3200 

Transverse cell pitch (S zo ) 

0.0147 m 

Input heat flux 

4, 8 and 16 KW / m 2 

Cell stut length (/) 

0.0147 m 

Inlet coolant temperature 

300 K 

Cell height (H) 

0.012 m 

Outlet coolant temperature 

305.0 K-360.0 K 

Material 

LM25 

Test section inlet pressure 

1 bar (ambient condition) 


Table of parameters 
for the current test 


Photograph of test rig 
and inserted LFM 
test sample 




































Governing Equations 


Momentum Equation 
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Energy Equation 
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Constant heat flux 

Flow and heat transfer are fully developed 
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= permeability 
p f = fluid density 
F = inertial coefficient 
s = porosity 

wetted area ratio 
k se , kj e = effective thermal 
conductivity 

Cf = specific heat of fluid 
h = interstitial heat transfer 
coefficient; 
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Results of Hydraulic Resistance 
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Infrared thermal image of top LFM substrate 
showing linear increase of temperature 
in constant heat flux boundary condition 
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Coolant velocity at exit 
Coolant velocity at inlet 
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Velocity (nVs) 


10 


Coolant velocity distributions along 
the channel height 

for both inlet and outlet of the test section 
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Comparisons 
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Comparison of overall Nusselt number with experimental data h- 

'» 1 fJ> 

UNIVERSITYOF- 

CAMBRIDGE 


Nu = ™*- 


'/ 



The Cambridge Centre for Micromechanics & Whittle Laboratory CUED 
























Interstitial heat transfer coefficient, h proposed by 
Zukauskas for Staggered banks of cylinder array 

Staggered cylinders 
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Effect of yaw angle 


h = 0.9 h 
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Re. =1-500 


= 500-1000 


= 10 3 -2xl0 : 
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Average Nusselt number 


Nu 
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Reynolds number 
based on 
cylinder bar 
diameter used, 
where 


h 


Nu d = 


k f / d 


Re. -U d/v 

a m 

and Prandtl number 
Pr is assumed to be 
constant e.g. 0.71 
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Comparisons with other heat sink media 

(from Kays and London, “Compact heat exchangers”) 



• 

LFM: Orientation A 


o 

LFM: Orientation B 


m - 

Bank of cylinders 


<0 

Wavy-fin plate-fin surface 
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Cross-rod matrices, inline stacking 



Plain pin-fin 


X 

Infinitely randomly stacked sphere matrix 
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Future Works 


(1) Detailed local heat transfer mechanism 
is anticipated by using Thermochromic 
Liquid Crystal (TLC) / Infra Red imagings. 


Heated airflow 

- > 



Liquid crystal coated 
Hollow cylinder 



Flat black coating 
Liquid Crystal 


Perspex 


Colour change of TLC is monitored 
from inside of the cylinder 
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(2) Construct heat sink design parameters 
using LFMs e.g. cell size effects, number 
of stacked layers, aspect ratio of the heat 
sink channels etc. 



Test section 


Traverse gears 



Test rig for single cylinder and a different 
version of rig for LFMs will be constructed. 
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Light weight metal foam heat exchangers 


-Examples for heat exchanger application 
-Flow resistance 


-Thermal resistance 
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Fabrication of metal foam heat exchangers 



Metal foam compact heat exchanger 
for high temperature service. 

Foam material is PFCT's FeCrAlY. 


SEM micrograph of a foam strut sintered 
to a solid tube. 

Bonding region shows metallurgical 
sintering between foam and solid. 



UNIVERSITY OF 

CAMBRIDGE The Cambridge Centre for Micromechanics & Whittle Laboratory 


CUED 








Fabrication of metal foam heat exchangers 



Example assemblies 
manufactured 
in a proprietary 
co-sintering technique 
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Governing Equations 

Momentum Equation 


0 = -V/> 


Energy Equation? 
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a = wetted area ratio 

k d = thermal dispersion conductivity 

= G>( Re * p hV M »A 

C D = 0.10 

Re t =(«„VZ)/r 

Pr«=l“C/)/*. 

k e = effective conductivity; kj e + k se 
F = inertial coefficient 

K = C(\-s) m (d f /d P ) n d 2 P 

where C= 0.00073, m = -0.224, n = -1.11 


Fully developed velocity field 
Constant heart flux q = q w = k fe 


h = 0.52Re 05 Pr 037 k f /d f 
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Results of Hydraulic Resistance 

- Static pressure drop per unit length - Friction Factor 



Operating range Re Dh = p UD h /p 


3000 < Re n <22000 
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a S-5: 5%, 60 PPI 



a S-6: 10%, 60 PPI 

• S-7: 7.5%, 30 PPI 
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Average Nusselt number. 


Case 1. fixed pore size as 30 PPI Case 2. fixed relative density as 10% 



AT =(YT 1/n-T h = q /A T and Nu = hD h /k f 

avg I iLj w,i I m A avg 
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Predictions of local and overall Nusselt numbers 
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Re = m h 
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Nu(x) = q Dh 


K~T m k f 





Re Re 


Variations of overall Nusselt numbers with Reynolds numbers for Porvair material 
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Nu(x ) = 


Z. 







X Re 


Boundary effects on local and overall Nusselt numbers 
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Comparisons with solid strut metal foams 

[from Calmidi et al. (2000)] 



5000 10000 15000 20000 25000 


Re 


Dh 


Thermal conductivity 

k = 16 W/mK (FeCrAlY) 

s 

~ 160 W/mK (Al alloy, 

T-6201) 


Nusselt number 

based on hydraulic diameter 

Nu = hD h /k f 

where red, green and yellow 
filled symbols indicate 
10%, 7.5% and 5% 
relative density, respectively 
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